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ABSTRACT 


The  Kalman  Filter  Is  Implemented  to  process  range  data  output  from 
the  Cubic  Model  40  Autotape  system,  a surface  position  locating  system 
currently  employed  on  the  underwater  tracking  ranges  at  Oabob  Bay  and 
Nanoose.  Results  are  presented  for  different  measurement  noise  and 
forcing  function  noise  statistics. 
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I.  INTRODUCTION 


The  Cubic  CM-40  Autotape  is  a microwave  distance  measuring  system 
used  (by  the  U.S.  Navy  at  its  acoustic  underwater  tracking  ranges  at 
Dabob  Bay  and  Nanoose)  to  provide  reference  position  information  for 
units  on  the  surface  and  in  the  air  above  the  range.  This  portable 
system  consists  basically  of  an  interrogator  which  is  operated  aboard 
the  unit  to  be  tracked,  two  responders  operated  at  two  different  shore 
sites  and  the  associated  antenna/RF  assemblies.  Required  support  systems 
l Include  a data  display  and  recording  setup  and  an  ADP  facility  for  off- 

line processing  of  the  Autotape  data.  Figure  1 shows  the  Autotape  system 
components  and  Figure  2 shows  a typical  application  geometry. 

Historically,  the  Autotape  has  been  used  in  such  applications  as 
tracking  hydrophone  array  survey,  buoy  and  hydrophone  array  planting  and 
as  a reference  position  Indicator  for  calibrating  other  position-finding 
devices  against.  Generally,  the  Autotape  has  been  used  where  an  extremely 
high  degree  of  accuracy  Is  not  required. 

In  operation,  the  system  will  provide  for  the  display  and  recording 
of  two  ranges  simultaneously,  once  per  second,  the  ranges  being  those 
between  the  Interrogator  and  each  of  the  responders.  The  ranges  are 
computed  from  the  phase  delay  between  the  output  of  the  modulation  signal 
generator  and  a signal  which  has  traveled  from  the  Interrogator  to  a 
responder  and  back.  Ranging  accuracy  Is  stated  by  the  manufacturer  to 
be  ± 0.5  meter  + 10  ppm  x range.  Ranging  frequencies  of  1500  KHZ,  150 
KHZ  and  165  KHZ  modulate  a 3000  MHZ  carrier,  yielding  a maximum  unambiguous 
range  of  10,000  meters  with  a resolution  of  0.1  meter.  However,  Independent 
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testing  by  the  U.S.  Navy  [Reference  1]  has  shown  that  system  accuracy  may 
not  be  quite  as  good  as  stated  by  the  manufacturer. 

The  accuracy  of  the  Autotape  system  Is  principally  dependent  upon 
range  errors,  the  geometry  of  the  system  and  the  method  of  data  reduction. 
These  factors  are.  In  turn,  affected  by  propagation  velocity,  system 
stability,  range  dependency,  land  survey  accuracy,  system  geometry, 
slope  reduction  and  data  smoothing.  A final  anomaly  which,  depending 
upon  the  application,  can  substantially  degrade  the  quality  of  the  data- 
stream  out  Is  the  orientation,  over  time,  of  the  Interrogator  antenna  in 
the  vertical  dimension.  The  Interrogator  antenna  has  only  a 10  degree 
vertical  beam  width.  Thus,  If  the  system  Is  being  used  on  a platform 
such  as  a moderately  maneuvering  helicopter  or  a ship  rolling  substan- 
tially In  the  seaway,  the  system  tends  to  frequently  lose  track,  resulting 
In  fairly  long  streams  of  useless  data. 

Present  data  reduction  techniques  employed  when  the  system  Is  used 
on  either  of  the  ranges  (Dabob  or  Nanoose)  employ  two  overall  Iterations. 
The  first,  or  initial  processing,  administers  the  following  three  correc- 
tions to  the  raw  range  data: 

1.  Range  Calibration  Correction:  This  is  a fixed  value  (meters)  added 
to  or  subtracted  from  each  range. 

2.  Propagation  Velocity  Correction:  This  Is  a variable  correction  due 
to  the  atmospheric  Index  of  refraction  at  the  particular  time  and 
place  of  the  exercise. 

3.  Slope  Reduction  Correction:  This  reduces  both  range  measurements 
(which  are  actually  slant  ranges  because  the  interrogator  and  the 
responders  are  not  normally  located  at  the  exact  same  elevation) 
to  a common  horizontal  plane  at  sea  level. 

Subsequent  processing  of  the  data  Includes  conversion  of  the  corrected 
ranges  to  a rectangular  x-y  range  coordinate  system  and  a moving  average 
smoothing  technique  which  employes  curve  fitting  algorithms  (linear. 
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parabolic  or  logarithmic)  to  reduce  the  data  to  Its  final  form.  Not 
uncommonly,  as  a result  of  the  total  reduction  effort,  the  net  remainder 
Is  an  Inadequate  data  package  (In  terms  of  quantity)  for  proper  final 
evaluation. 

Figure  3 Is  a rectangular  plot  of  the  raw  ranges  recorded  during  a 
recent  array  survey.  The  purpose  of  this  project  has  been  to  design  a 
filter,  a Kalman  filter,  which  would  provide  more  accurate  range  data, 
as  well  as  one  that  would  track  through  the  periods  of  "lost  track" 
ranging,  thereby  providing  a significantly  larger  final  volume  of  data 
for  evaluation.  This  paper  presents  the  basic  theory  necessary  and 
Includes  the  final  version  of  the  filter. 
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II.  THE  FILTER  THEORY  AND  DESIGN 


A,.  THE  SYSTEM  DYNAMIC  MODEL 

A common  application  for  the  Autotape  system  is  Its  use  as  a reference 
position  locator  on  the  surface  unit  conducting  an  acoustic  hydrophone 
array  (range)  survey.  The  usual  exercise  plan  will  call  for  a service 
unit,  carrying  the  interrogator  and  equipped  with  an  acoustic  plnger 
mounted  on  the  underwater  hull,  to  transit  three  concentric  circular 
tracks,  centered  above  the  array,  with  track  radii  ranging  from  100  to 
1,000  meters,  at  speeds  of  up  to  eight  knots.  The  direction  of  rotation 
for  the  outer  track  will  normally  be  opposite  to  that  of  the  middle 
circle.  While  the  service  unit  is  being  tracked  via  Autotape,  It  Is  also 
being  tracked  by  the  acoustic  array.  By  comparing  the  acoustic  position 
data  with  that  from  the  Autotape,  a digital  computer  is  able  to  compute 
actual  position  and  attitude  of  the  array. 

The  desired  estimates  will  be  those  of  position  and  velocity,  Rl,  R2, 

• • 

Rl,  R2.  It  is  proper  at  this  point  to  define  a number  of  terms  and  to 
summarize  some  pertinent  results  of  observer  theory.  First,  we  may  define 
a fourth  order  state  vector: 


Recall  that  a linear  system  can  be  described  in  the  continuous  time 


domain  as: 


* (t)  * Ax  (t)  + D w (t) 


where: 


x(t)  Is  the  n -element  column  vector  of  the  states 

A and  0 are  nxn  and  nxp  matrices  describing  system  dynamics 

w(t)  Is  a q-element  vector  of  random  noise  Inputs  to  the 
system 

The  system  measurements  may  be  expressed  as: 
z(t)  x (t)  + v^  (t) 

where:  z.(t)  Is  the  q-element  vector  of  system  measurements 

Is  the  qxn  weighting  matrix  for  the  measurements 
v(t)  Is  the  q-element  vector  of  random  measurement  noise 

The  corresponding  linear  discrete  model  may  be  written  as: 
x(k  + 1)  * g.x.(k)  + r w(k) 

with  no  deterministic  Inputs  to  the  system. 

Also,  ,z(k)  3 H x(k)  + v(k) 


For  the  system  under  consideration.  It  can  be  shown  that  the  state 
transition  matrix 


£ 


10  10 
0 10  1 
0 0 10 


for  a sampling  Interval  T of  1 second.  A block  diagram  of  the  system  Is 
shown  In  Figure  4. 
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The  following  assumptions  will  be  made  regarding  the  noise  processes 
and  the  initial  state,  xjto)  of  the  plant  [Ref.  2]: 


The  measurement  noise  has  zero  mean.  Is  uncorrelated,  and 
E [ v.(k)  _vT( j ) 3 ■ R(k)  Jy,  where  % Is  the  kronecker  delta 

The  forcing  noise  has  zero  mean.  Is  uncorrelated,  and 
E [ w(k)  wT(j)J  * £(k)  *kj 

The  forcing  noise  and  measurement  noise  are  uncorrelated. 

The  Initial  state  Is  a random  variable  with  known  mean  and  covariance, 
and 

E C f x(o)  - x0}  fx(o)  - x0)  T]  » PQ 

The  measurement  noise  and  Initial  state  are  uncorrelated. 

The  forcing  noise  and  initial  state  are  uncorrelated. 

The  Kalman  Filter  equations  and  their  derivation  are  well  known 


[Ref.  2],  [Ref.  31: 

6(k)  » P(k/k-l)  HT(k)  [H(k)  P(k/k-l)  HT(k)  + ^(k)]"1  (1) 

P(k/k-l)  - £ P(k-l/k-l)  £T  + £ (2) 

P(k/k)  - H - G(k)  H(k)3  P(k/k-l)  (3) 

£(k/k)  - x(k/k-l)  + G(k)  [z(k)  - H(k)  x(k/k-l)]  (4) 

|(k/k-l)  » £(k/k-l)  x(k-l/k-l)  + r(k/k-l)  w(k-l)  (5) 

% 


Where  the  notation  (k/k-1)  Interprets  as  the  value  of  the  parameter  of 
note  at  time  k given  measurements  at  times  up  to  and  Including  time  k-1. 
(k/k)  and  (k-l/k-1)  have  similar  Interpretations.  The  % denotes  the 
estimate  of  x. 
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G(k)  represents  the  filter  gain  at  time  k.  P represents  the  covariance 
of  estimation  error; 


T ) e?(k/k) 

P (k/k)  » E [e(k/k)  eT(k/k)]  - E / ,z  [e^k/k)  e2(k/k)~ 


e^k/k) 


en(k/k) 


-en(k/k)] 


ej(k/k) 


e1(k/k)  e2(k/k)  — e^k/k)  en(k/k) 


e2(k/k)  e1(k/k)  e|(k/k) 


— e2(k/k)  en(k/k) 


en(k/k)  e1(k/k)  en(k/k)  e2(k/k)  — eZ(k/k) 


where  e(k/k)  * $(k/k)  - x(k).  A complete  standard  block  diagram  for  the 
filter  and  an  Information  flow  diagram  are  Included  as  Figures  5 and  6 
as  slightly  different  viewpoints  from  which  the  system  may  be  viewed  and 
understood.  Figure  7 shows  a timing  diagram  of  the  various  quantities 
contained  In  the  filter  equations. 
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H(k-l)  R(k-l) 


£(k-l)  £(k-l) 


Compute 

| 

Compute 

P(k/k)  - [I-6(k)H(k)]  P(k/k-l) 

■ 

P(k/k-l)  - l P(k-l/k-l)JiT+  ^ 

H(k-1),  R(k-1) 

A 


$(k-l/R-2) 


P(k-l/k-2) 


* (k-l/k-1) 


H(k),  R(k) 


$ (k/k-1) 


I(k-l),  5(k-l) 


P (k-l/k-1) 


P (k/k-1) 


x (k/k) 


£(k),  £(k) 

» 


P (k/k) 


Lk-1 


time 


FIGURE  7:  Timing  Diagram  of  Filter  Equation  Quantities 
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B.  THE  PROCESSOR 


Appendix  A Is  a flowchart  of  the  Kalman  filter  program  utilized. 
Initially,  the  matrices  describing  the  physical  system,  the  noise 
statistics  and  other  program  parameters  are  read  Into  storage  and  printed 
out.  The  discrete  state-transition  matrix.  Phi,  Is  computed  and  printed 
out  and  the  gain  schedule  Is  computed  and  printed  out.  It  Is  seen  that 
the  elements  of  the  gain  matrix  reach  a steady  state,  and,  for  example, 
with  both  the  R and  £ matrices  being  Identity  matrices,  the  gain  reaches 
steady  state  between  k*5  and  k«10.  Therefore,  In  the  main  Iteration  loop, 
the  filter  will  essentially  be  a constant  gain  filter  for  k>10. 

Next,  the  main  Iteration  loop  commences.  The  Initial  measurements 
are  read  and  xl(0/-l)  and  x2(0/-l)  are  Initialized  to  these  values. 
x3(0/-l)  and  x4(0/-l),  representing  the  rates,  are  set  to  the  mean 
constant  value  (In  the  respective  directions)  of  4.0  meters  per  second. 

The  Autotape  output  Is  a 5 significant  figure  output,  modulo  10,000, 
reading  to  0.1  meter.  Inherent  in  the  output  is  a major  degree  of  jitter 
In  the  two  most  significant  digits,  which  would  significantly  distort 
the  covariance  of  measurement  noise.  Therefore,  as  an  option,  measure- 
ments could  be  gated,  and  the  gain  automatically  set  to  zero  In  those 
cases  where  the  residue  falls  outside  of  a maximum  reasonable  bound. 

Commencing  with  k«0,  and  utilizing  the  known  values  for  x(0/-l)  and 
P(0/-1),  the  Kalman  filter  equations  are  solved  Iteratively  In  the 
following  manner  [see  page  15,  equations  (l)-(5)]: 

(1),  (3),  (4), 

Increment  k to  k-1 
(5),  (2),  (1),  (3),  (4), 

Increment  k to  k«2 
(5),  (2),  (1),  (3),  (4), 
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Also  computed  on  each  Iteration  are  the  error  residues: 

RES  * z - x(k/k-l) 
and  the  one-step  prediction  errors: 

ERR  = x(k/k)  - x(k/k-I) 

Finally,  the  computations  are  tabulated  and  plots  are  produced. 

C.  NOISE  AND  ERROR  CONSIDERATIONS 

Reference  1 documents  an  Autotape  evaluation  which  was  conducted  In 
1971.  The  error  geometry  Is  shown  In  Figure  8.  Graphically,  position 
Is  determined  by  locating  the  crossing  point  of  the  two  range  arcs.  In 
conjunction  with  a knowledge  of  the  baseline  formed  by  the  two  responders. 
Since  each  range  has  an  associated  standard  deviation  (error),  the  point 
can  actually  be  enclosed  In  a parallelogram  which  defines  the  probable 
position  within  one  standard  deviation  of  the  ranges.  The  shape  of  the 
parallelogram  will  vary  with  the  position  of  the  crossing  point  relative 
to  the  baseline,  as  Indicated  In  Figure  8.  It  can  be  shown  that  the 
maximum  probable  error  (MPE)  will  be  minimized  where  the  range  arcs  are 
orthogonal.  Figure  9 diagrams  error  contours  which  are  actually  the 
locll  of  constant  WE  for  two  particular  responder  sites  on  the  Nanoose 
Range.  Table  1 summarizes  pertinent  results  of  the  study. 
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resent  maximum  probable  positional  error  in  feet.  End  points  of  arcs 


TABLE  1 


Average  Range  Errors  (feet) 


R-l 

R-2 

PI 

No.  Points 

Error 

Average 

Standard 

Deviation 

Error 

Average 

Standard 

Deviation 

Array  04 

30 

- 0.5 

2.8 

m ■ 

2.8 

Array  07 

49 

- 1.3 

2.3 

' ■ BiW 

2.4 

Array  08 

10 

- 0.8 

4.4 

2.8 

Array  09 

25 

3.8 

2.6 

i 

2.2 

Average 

0.3 

3.0 

- 

2.6 

For  the  purpose  of  modeling  the  covariance  of  excitation  noise.  It 
was  assumed  that  the  service  unit  transited  an  800  meter  circle  at  an 
average  speed  of  eight  knots.  Then: 


Filter  performance  was  Investigated  for  2.  = I»  and  -01I» 


where  the  a priori  )<(0/-l)  Is  known  to  be  a reasonably  good  estimate  — 
approximately  the  same  accuracy  as  an  observation. 


D.  PROCESSOR  PERFORMANCE;  AUTHOR'S  CONCLUSIONS 

Table  2 summarizes  a comparison  of  the  Kalman  filter  performance  with 
the  results  of  the  (corrected)  processing  by  the  program  presently  being 
used  for  the  cases  2 * I»  “ i*  and  2 * 0.1I_,  R s I_.  Figures  10,  11, 

12  and  13  are  residue  and  error  plots  for  the  example  2 * .01I_,  R * K 
It  Is  seen  that  the  Kalman  filter  will  satisfactorily  handle  the  data 
where  the  measurement  noise  statistics  approximate  those  used  In  the 
model.  However,  for  the  noise  resulting  from  the  jitter  which  appears 
In  the  "hundreds"  and  "thousands"  digits,  the  filter,  as  configured 
without  a gate,  will  estimate  with  considerable  error.  The  raw  range 
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R2  was  clean  of  this  particular  noise  element,  and  the  results  as 
Indicated  by  Figures  12  and  13  were  superior  to  those  for  Rl. 

It  Is  suggested  that  the  Kalman  filter  be  used  as  the  first  Iteration 
processing  of  the  Autotape  output. 
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TABLE  2 

TABULATED  PROCESSOR  COMPARISON 


RESIDUE  t <X10‘ 


ERROR  1 (XIO*) 


III.  FUTURE  FILTER  IMPROVEMENTS 


The  filter,  as  designed,  will  process  by  off-line  (forward)  filtering 
of  the  range  measurements.  It  Is  suggested  that,  as  an  effort  to  further 
Improve  upon  the  quality  of  the  processed  data,  a fixed-interval  smoothing 
algorithm  (the  initial  and  final  times,  0 and  T,  are  fixed,  and  the 
estimate  ^ (t/T)  is  sought)  be  Incorporated. 

For  the  system  and  measurements  described  by: 

*.  * Fx  + Gw 
z * Hx  + v 

the  equations  defining  the  forward  filter  are,  in  the  time  domain  [Ref. 3]: 


i * f*  + PhV1  [z  - Hx],  x « ^ (1) 

P = _FP  + PFT  + GQG1  - PHYV  P(0)  - ^ (2) 

To  write  the  backward  filter  equations,  setT»  T - t.  Then  -gjp-  ■ -g|-. 


Also, 


-g=r  * -Fx  -Gw,  for  0 < T <;  T,  denoting  differentiation  with 
~ respect  to  backward  time. 


,z(T)  * Hx  + v,. 


Then,  by  analogy,  the  backward  filter  equations  can  be  written  by 
changing  £ to  -£  and  £ to  -£,  resulting  In: 

-jri,  ■ -fs,  ♦ f/i'1  tz  - jgJ,] 

•no  -pb  • -ffb  - -VT  * ■ -V-VX  (3 


forward  filter 


i 

\ 

\ 

i 


From  Figure  14,  It  can  be  seen  that  the  smoothed  estimate  at  t1me*T 
must  be  the  same  as  the  forward  filter  estimate  at  that  point,  l.e., 

x (T/T)  - x (T) 
and  P (T/T)  - P (T) 

which  yields  the  required  boundary  condition  on  PjJ1, 

Pjj1  (t-T)  » 0,  or  P^1  (TfcO)  - 0 (4) 

with  the  boundary  condition  on  ^ (T)  not  yet  known.  Therefore,  define 
the  new  variable: 


I 


l 


f 


f 

i 

f 

I 


Reformulation  In  terms  of  P^  yields: 

= -51  41 

Thus,  equation  (3)  can  be  written  as: 


-J r 51  ■ f;1  £ + ft  g - p^1  a asT  p^1  + HT  R'1  H (7) 


for  which  equation  (4)  Is  the  appropriate  boundary  condition. 


Differentiating  equation  (5)  with  respect  to  *T,  and  with  some 
substitution  and  manipulation,  we  arrive  at: 

-jjp-  s * | FT  - Pjj1  G a iT  J s + HT  R"1  z (8) 

for  which  equation  (6)  Is  the  appropriate  boundary  condition.  Equations 
(1).  (2),  (7)  and  (8),  along  with: 


P'1  (t/T)  - P"1  (t)  + P^1  (t) 
x (t/T)  - P (t/T)  [P"1  (t)  x (t)  + P^1  (t)  (t)] 

define  the  optimal  smoother. 


Many  forms  of  the  smoothing  equations  may  be  derived.  The  form 
proposed  for  use  In  this  particular  case  Is  the  Rauch-Tung-Strlebel  form, 
with  the  discrete-time  expressions  summarized  as  follows: 


i 


Smoothed  State  Estimate  £(k/N)  *£(k/k)  + [£(k+l/N)  -£(k+l/k)] 


r 


where 

^ = P(k/k)  &(k)T  P(k+l/k)-1 
for  k * N-l 


Error  Covariance 
Matrix  Propagation 


P(k/N)  - P(k/k)  + [P(k+1/N)  - P(k+l/k)]  aJ 

also  for  k a N-l 


Solution  of  the  equations  would  proceed  as  follows;  As  an  example, 
and  because  it  Is  slightly  easier  to  see  when  actual  times  are  used, 
suppose  NN  * 100.  On  the  forward  filter  pass,  the  values  of  £(k/k), 
£(k/k-l),  P(k/k)  and  P.(k/k-l)  would  be  computed  and  stored.  On  the 
final  iteration  of  the  forward  pass,  with  K * NN  * 100, 

£(100/100)  *£(100/99)  + G(100)  [z(100)  - Hxt 100/99)] 
l.e.,  we  have  computed  and  stored  £(100/100). 

Now,  the  smoothing  process  comnences  In  the  reverse  direction. 
Decrement  k to  k * NN-1  * 99,  then 


£(99/100)  *2(99/99)  + A(99)  [£(100/100)  - £(  100/99 )] 


stored  stored  stored 


and  A (99)  - P(99/99)  £T  P(  100/99  f1 

stored  stored 
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let  k - NN-2  - 98,  then 


r 


"5(98/100)  « 5(98/98)  + A(98)  [5(99/100)  - x( 99/98)] 

stored  computed  stored 

last 

Iteration 


and  A(98)  - £(98/98)  £T  P(99/98)-1 
stored  stored 


Also,  for  each  of  the  two  preceding  iterations. 


£(99/100)  * £(99/99)  + A(99)  [£(100/100)  - £(100/99)]  AT(99) 
stored  computed  stored  stored 


P_(98/100)  * P( 98/98)  + A(98)  [£(99/100)  - £(99/98))  AT  (98) 
stored  computed  computed  stored  computed 


etc. 


It  Is  seen  that  the  smoothing  process  does  not  involve  the  processing 
of  actual  measurement  data.  It  does,  however,  utilize  the  complete 
filtering  solution,  and  so  fixed  interval  smoothing  cannot  be  done  real- 
time, on-line.  It  must  be  done  after  all  the  measurement  data  are 
collected.  Consequently,  computation  speed  will  not  be  the  most  Important 
factor.  Storage  requirements  could,  however,  conceivably  be.  In  that  the 
quantities  to  be  stored  on  the  forward  pass  are  arrays.  It  is  seen  that, 
should  an  exercise  run  in  excess  of  30  minutes,  retention  of  the  data  at 
each  mark  could  require  In  excess  of  100K  bytes  of  memory,  which  could 
limit  the  facilities  upon  which  the  processor  could  be  utilized. 
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A 


Call  PROD 
xCOR  = G zMOD 
= G(z-H  x(k/k-l) 


COMPUTER  OUTPUT 

i “ I*  B. » i 


3AM  32 

4232. 2 

4973.1 
4372*7 

4372.4 

4363.8 
4 352  « 2 

4952.7 
4353. G 

4252.5 

4951.7 
4=49.7 

4342.7 

4242.3 
4 342 .5 

4337.2 

4333.2 
4931.  e 

4323.4 

4922.5 

4921.1 

4312.5 

4215.2 
422?.? 

4312. 2 

4224.6 

4322.7 

4623.2 

4336.7 

4622. 7 

4337.2 

4332.2 
4332  .2 

4361.2 

4377.7 

4271.9 
4 372.2 
4f£?  .? 
43c". 4 

4322.6 
4 3 i 3 . * 

4827.2 
4351*2 
4E43.6 

4347.7 

4344.1 

4359.3 

4341.6 

4452.1 

4322 .2 

4325.9 
4S24.C 


I1.7EP.E2  7 1 

FUTEPZ3  F 

4622.4 

4232.2 

4627.4 

433C.4 

4U  w^. 

4972.2 

4 52  5.3 

4972*2 

4612.6 

4363. C 

4S4C.5 

4962*9 

4646  .C 

4962*4 

4&49. 1 

4359.7 

5474 .2 

4955.2 

4373.1 

4321 *2 

46  81.8 

4942.4 

4622.4 

4944*2 

4 C 4 2 .3 

424C.7 

4552*7 

4341*3 

4662.7 

4927.8 

4S73.7 

4923*4 

4614.6 

4931*2 

4621.2 

4928*4 

4622.9 

4923* C 

4673.5 

4322  .7 

4614.4 

491S-2 

4652.3 

4316. 2 

47C2.7 

4912*5 

4714.4 

4923*2 

4712.5 

49C5  *2 

4712.3 

4 322  *5 

4716.2 

4327*7 

4722.7 

4335*3 

4716.2 

4322.2 

4723  .7 

4333*2 

4732.C 

4662*1 

47*6  .7 

4332*2 

4 7“2  .2 

4826.2 

4 74  5 *6 

4377*3 

4745.2 

4872.5 

4772.5 

4 372  .3 

4727.9 

4862.1 

4762.5 

4353.7 

4 7C  4 .2 

4362*2 

.4777.5 

4373.1 

4772  .2 

4355.6 

4773.7 

4352*2 

4714.2 

4343.6 

4 / C2  • 2 

4346*9 

4721.3 

4344.2 

4723.5 

4J4C  *1 

4714.3 

4842.7 

4713.2 

4333.4 

4787  .2 

4329.2 

4323.2 

4326.5 

4312. 2 

4323*8 
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£ - i,  r - 1 


RES I3UE 

•0 


RESIDUE 

•C 


10.9 

-7.1 

s.c 

-9.1 

1.2 

1.3 

.7 

-1.1 

-3.2 

-4.1 

4.1 

4.7 

-l.C 

-.3 

595.6 

-1.5 

# **  •* 

-l.C 

iie.c 

1.4 

124.9 

-2.7 

127.4 

• 3 

72.3 

5.3 

25*6 

-2.9 

3.3 

-2.2 

-76.6 

1.7 

-7E.S 
13  .5 
2E.C 
19.8 

e.s 

1.7 

61.4 

-20.7 


•* 
2*2 
-.3 
-1.7 
4.9 
-7.3 
• 2 
2*1 
• 9 


ERROR 

• C 

9.0 

4.1 

1.0 

.6 

-2.6 

3*3 

-.8 

£21*6 

-37.  C 
1C2«7 
204-7 
59.5 
21. C 
6.3 
-63.0 


-E2.2 

15.2 
23. C 

16.2 
7.C 
1.4 

86.9 

-17.0 


ERROR  2 

TIME 

•C 

1C5538 

-5.8 

105539 

-7.5 

1C554C 

1.1 

1C5341 

-.9 

IC554Z 

-3.3 

105543 

2.9 

1C5544 

-.3 

1C554S 

-1.2 

. 1C5546 

-.3 

105547 

1.2 

1C5S48 

-2.3 

105349 

.6 

1C555C 

4.7 

ICS351 

-2.4 

105552 

-1.3 

105553 

1.4 


2C  *0 

-.  1 

16.4 

-.1 

2 6*5 

-3.1 

23.4 

-2.5 

13  .6 

2.2 

15.3 

1.3 

$-4 

2.1 

7.7 

1 .7 

35.4 

-.9 

70-2 

-.7 

-22.3 

-3.7 

-18.3 

-3.C 

-25.7 

3.3 

-21.1 

3.1 

-20. £ 

-2-1 

-16.9 

— 1*7 

-13.2 

-1.1 

-1C.  8 

-.9 

.4 

1.5 

.2 

1.2 

2*2 

2.3 

1.8 

1.3 

-2.6 

-1.2 

-Z.2 

-l.C 

-1.3 

-2.4 

-1.1 

-2  • C 

1.4 

-1.2 

1.1 

-l.c 

• 1 

4.3 

.1 

3.5 

2.C 

1.6 

1.7 

1*3 

• C 

-1.5 

• C 

-1.2 

-2.7 

-3.3 

-2.2 

-2  *7 

1.3 

2.4 

1*1 

1.9 

.5 

2.C 

.4 

1*6 

*6 

-1*6 

. 5 

-1.3 

-2.5 

-3.3 

-2.9 

-2.7 

.7 

2.5 

• 6 

2.1 

1.8 

2.1 

1.5 

1*7 

-3.1 

-2.9 

-2.5 

-2.4 

1 .8 
-•2 
'1  .4 
4.C 
'5*0 
•2 
1.7 

.e 


IC5554 

1CS555 

1C55E6 

105537 

105558 

1C5559 

1C5GCC 

105521 

1C5602 

105623 

1C5GC4 

1C5325 

1C5GC6 

1C3S27 

1C56C8 

1C5S29 

1C5G1C 

1C5S11 

1CSG1Z 

103613 

105614 

105515 

1C561C 

105317 

105618 

105319 

1C5C2C 

1C3321 

1C5622 

105623 

105624 

105625 

1C5626 

1C5327 

105628 
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2 « o.n,  r * i 


K 

RAW  Rl 

RAW  R2 

FILTERED  Rl 

FILTERED 

C 

4622*4 

4982.2 

4622*4 

4982*2 

1 

4629*3 

4979.1 

4624*7 

4982*1 

2 

4633*C 

4972*3 

4628*8 

4977.5 

3 

4636*0 

4972.4 

4633.1 

4974*6 

4 

4639*7 

4968*8 

4637*6 

4970*8 

5 

4639*9 

4962.2 

4640*1 

4965*0 

E 

4646*7 

4963.2 

4645.  C 

4962*5 

7 

4648*9 

4959.6 

4648*7 

4959*6 

8 

5652*3 

4955*6 

5230.4 

4956.0 

9 

4656*6 

4351.7 

4986*8 

4952*2 

1C 

4660*8 

4949.7 

4818*7 

4949.3 

11 

4657.7 

4943.7 

4713.5 

4944.7 

12 

4666*5 

494C.8 

4662*8 

4940. 9 

13 

4673*6 

4942*6 

4646.2 

4940.3 

14 

4674*3 

4937.3 

4645.3 

4937.5 

15 

4676*  8 

4933. C 

4652.3 

4933*7 

IE 

4601*0 

4931.E 

4616*5 

4931*2 

17 

4604*6 

4928.4 

4600* 2 

4928*4 

18 

4607* 9 

4922*5 

4596. 2 

4923*8 

19 

4611*3 

4921*1 

4599.2 

4920*8 

2C 

4616-1 

4919*6 

4605.5 

4918*7 

21 

4702-0 

4916. C 

4660.0 

4915*9 

22 

4702*7 

4909.8 

4692.0 

4911.2 

23 

4709*8 

4910.0 

4711.8 

4909. 0 

24 

4711-8 

4904.6 

4720* 9 

4905*2 

25 

4711*6 

4900.3 

4722.8 

4901.0 

26 

4718*1 

4898.0 

4725.0 

4897.7 

27 

4724*1 

4898.7 

4728.0 

4895.7 

28 

4725*7 

4892.7 

4729.4 

4892.7 

29 

4728* S 

4887.6 

4730.8 

4888.5 

30 

4733*2 

4882.9 

4733.7 

4883.8 

31 

4736*7 

4883.0 

4736.8 

4881.7 

32 

4742*4 

4881.2 

4741.4 

4880.1 

33 

4746*6 

4877.5 

4746.0 

4877.4 

34 

4748*5 

4871*9 

4749.1 

4873.1 

35 

4753*7 

4870.9 

4753.3 

4870.5 

3E 

4757*9 

4869.5 

4757.6 

4868.6 

37 

4762*6 

4865.4 

4762.2 

4865-7 

38 

4763*4 

4859.6 

4764.7 

4861.0 

39 

4767*7 

4858.5 

4767.9 

4858.1 

4C 

4772-6 

4857*0 

4772.1 

4856.1 

41 

4773*2 

4851.5 

4774.3 

4852 *3 

42 

4700*8 

4848*6 

4733.2 

4848.8 

43 

4703*  6 

4847.3 

4710.8 

4846.6 

44 

4706*9 

4844.1 

4701.8 

4843.9 

45 

4712*1 

4339.8 

4702. 0 

4840.3 

46 

4715*8 

4841.6 

4706.3 

4839.8 

47 

4718*3 

4332*1 

4711.6 

4834*4 

48 

4802*4 

4829.2 

4763.8 

4820*0 

49 

4804*8 

4826*9 

4795*0 

4826*6 

5C 

4807*8 

4824*0 

4812*0 

4823.6 

47 


4-o,n.  R-i 


RESIDUE  1 

RESIDUE  2 

ERROR  1 

ERROR  2 

TIME 

• C 

•0 

-0 

•C 

1C5538 

10  *9 

-7.1 

6.3 

-4.1 

105539 

1C-1 

-12.3 

5.6 

-7.1 

105540 

6.9 

-5.1 

4.0 

-3.0 

105541 

4.  9 

-4.7 

2.8 

-2.7 

1C5542 

-.5 

-6.6 

-.3 

-3.8 

1C5S43 

4-C 

1.6 

2*3 

• 9 

1C5S44 

.4 

• 1 

• 2 

• 1 

105545 

1CC0.0 

-.9 

578.2 

-.5 

105546 

-792*8 

-1.1 

-452.5 

— «6 

105547 

-374.2 

1.0 

-216.3 

.6 

1C5548 

-132.3 

-2.4 

-76.5 

-1.4 

105549 

e.s 

-.2 

5.1 

-.1 

10555C 

64.9 

S.5 

37.5 

3.2 

1C5551 

66.8 

-.4 

39-8 

-.2 

1C5552 

58. 1 

-1.8 

33*6 

-1.0 

105553 

-36.6 

.9 

-21.2 

.5 

105554 

10.3 

• 1 

6 *C 

.0 

1C5SS5 

27.7 

-3.0 

16.  C 

-1.7 

105556 

30.0 

.8 

17-3 

• 5 

1C5S57 

25.1 

2.2 

14.5 

1.2 

1C5558 

99.6 

• 2 

57.6 

• 1 

1CS559 

25.3 

-3.3 

14.6 

-1.9 

1C56CC 

-4.8 

2.3 

-2.8 

1.3 

105621 

-21.6 

-1.4 

-12.5 

-.8 

105602 

-26.5 

-1.6 

-15.3 

-.9 

105623 

-16.4 

.7 

-9.5 

.4 

1056C4 

-9.3 

2.5 

-5.4 

1.4 

105625 

-6.8 

.0 

-5*1 

• 0 

105606 

-5.5 

-2.1 

-3-2 

-1.2 

105627 

-1.1 

-2.2 

-.7 

-1.3 

1C56C8 

-.3 

3.0 

-•2 

1.8 

1CS629 

2-4 

2.7 

1.4 

1.6 

1C  56 1C 

1.5 

• 1 

.9 

.1 

1CSS11 

-1.5 

-2.9 

-.9 

-1.7 

1C5612 

.9 

1.0 

.5 

*3 

105613 

.7 

2.1 

• 4 

1.2 

105614 

1.0 

"•  6 

.6 

-.4 

1C561S 

-3.0 

-3.3 

-1.7 

-1.9 

105616 

* *6 

.9 

-•3 

.5 

105617 

1*2 

2.1 

.7 

1.2 

105618 

-2  .6 

-1.8 

-1.5 

-1.0 

105619 

-76.7 

-.5 

-44.3 

-.3 

10562C 

-17.0 

1.8 

-9.8 

1.0 

ICS621 

12*2 

.5 

7.C 

• 3 

1C5622 

23*9 

-1.3 

13.8 

-.7 

1CS623 

22.4 

4.3 

13.  C 

2.5 

105624 

16.0 

-5.5 

9*2 

“3*2 

105325 

91.6 

-1.9 

53. C 

-l.i 

1C5626 

22*8 

• 6 

13*2 

•4 

1C5627 

-5.9 

• 9 

-5.7 

• 5 

j 

ICS628 
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r 


2 * o.on,  r - i 


K 

RAW  Rl 

RAW  R2 

FILTERED  Rl 

FILTERED  R2 

0 

4622*4 

4982*2 

4622.4 

4932.2 

1 

4629*3 

4979.1 

4622*4 

4983*6 

2 

4633*0 

4972*3 

4624.3 

4981.6 

3 

4636*0 

4972*4 

4627.4 

4979.6 

4 

4639*7 

4968*8 

4631.3 

4976*5 

k 

4639*9 

4962*2 

4634.5 

4971*5 

6 

4646*7 

4963.2 

4639*5 

4967.9 

7 

4648*9 

4959*6 

4644.0 

4964.0 

8 

5652*  3 

4955*6 

5017.2 

4959.7 

9 

4656*6 

4951*7 

4936.2 

4955*2 

1C 

4660*8 

4949.7 

4864.3 

4951*4 

11 

4657*7 

4943.7 

4801*7 

4946.6 

12 

4666*5 

434C.8 

4763.9 

4942-3 

13 

4673*6 

4942*6 

4719.4 

4940*1 

14 

4674*3 

4937.3 

4694.3 

4937.0 

15 

4676*8 

4933.C 

4677.8 

4933*4 

16 

4601*0 

4931*6 

4639.3 

4930. 7 

17 

46C4.6 

4928*4 

4613*3 

4927.8 

13 

4607*9 

4922*5 

4597.4 

4923.9 

IS 

4611*8 

4921.1 

4589.6 

4920.8 

20 

4616*1 

4919.6 

4588.1 

4918.3 

21 

4702*0 

4916 .0 

4621.0 

4915.5 

22 

4702*7 

4909*8 

4648.5 

4911.5 

23 

47C9.8 

4910. C 

4672.7 

4908* 9 

24 

4711*8 

4904.6 

4691.8 

4905.3 

25 

4711*6 

4900*3 

4705.3 

49C1.5 

26 

4718*1 

4898.0 

4716.7 

4898.1 

27 

4724*1 

4886*7 

4726.2 

4895.5 

28 

4725*7 

4892.7 

4732.7 

4892*4 

29 

4728*5 

4887*6 

4737.2 

4888.6 

3C 

4733*2 

4382*9 

4741.1 

4884.5 

31 

4736*7 

4883.0 

4744.3 

4881.7 

32 

4742*4 

4381.2 

4747.7 

4879.4 
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